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ABSTRACT 


A  control  system  based  on  sliding  mode  control  and  the  linear  quadratic 
regulator  is  designed  to  stabilize  the  straight  line  motions  of  towed  vessels.  The 
control  technique  is  through  athwartship  movement  of  the  towline  attachment 
point  on  the  towed  vessel.  Control  design  is  based  on  the  linearized  sway  and 
yaw  equations  of  motion.  Numerical  simulations  for  both  the  linearized  and  the 
nonlinear  system  are  performed  and  demonstrate  the  added  robustness  of  the 
control  technique  employed. 
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I.  INTRODUCTION 


The  horizontal  plane  stability  of  towed  vessels  is  a  very  important  field.  If 
the  motion  of  a  vessel  imder  tow  is  imstable,  a  very  dangerous  situation  exists; 
the  motion  of  a  several  thousand  ton  vessel  could  endanger  the  ship  towing  it 
[Refs.  1, 2,  3].  A  method  needs  to  be  developed  not  only  to  ensure  the  stability  of 
the  towed  vessel's  motion,  but  to  control  it  [Refs.  4, 5].  If  the  motion  of  the  towed 
vessel  can  be  controlled,  the  next  step  is  to  optinuze  its  control,  both  for  safety 
and  fuel  considerations.  Usually  the  only  means  of  control  of  a  vessel  under  tow 
is  through  the  towline.  Under  these  conditions,  stability  of  the  vessel  is  marginal; 
in  some  cases,  the  vessel  is  actually  unstable. 

This  thesis  will  develop  the  equations  of  motion  for  a  vessel  in  tow.  It  will 
also  state  and  describe  the  unique  forces  associated  with  a  vessel  being  towed. 
Finally,  it  will  describe  how  state-space  theory  is  being  applied  to  control  the 
towed  vessel's  motion  through  the  movement  of  the  attachment  point  of  the 
towline  along  the  beam. 

A  control  system  based  on  sliding  mode  control  and  the  linear  quadratic 
regulator  for  the  linearized  equations  of  motion  will  be  presented,  with  results  in 
non-dimensional  form.  The  same  control  system  for  the  non-linearized  equations 
of  motion  will  also  be  presented,  along  with  the  results.  Comparison  of  results  for 
the  linearized  and  non-linearized  equations  of  motion  will  help  determine  the 
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robustness  of  the  system.  Matrix^  will  be  used  to  calculate  the  required  gains  for 
the  system,  and  the  simulation  will  be  done  using  a  Fortran  program. 
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II.  EQUATIONS  OF  MOTION 


The  two  equations  of  motion  will  be  developed  using  a  body-fixed  frame  of 
reference.  The  origin  of  this  reference  frame  is  at  the  center  of  gravity  of  the 
towed  vessel.  In  general,  the  vessel  has  six  degrees  of  freedom: 

1.  surge  along  the  Xg-axis 

2.  roll  about  the  Xg-axis 

3.  sway  along  the  yg-axis 

4.  pitching  about  the  yg-axis 

5.  heave  along  the  Zg-axis 

6.  yaw  about  the  Zg-axis 

This  thesis  will  only  deal  with  motion  in  the  horizontal  plane.  Also,  a  constant 
surge  velocity  along  the  Xg-axis  will  be  assumed.  Other  assumptions  include: 

1.  no  towed  vessel /towing  vessel  interaction 

2.  massless,  inextensible  towline 

3.  no  wind,  wave,  or  current  disturbances 

4.  the  body  is  symmetric 
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Since  motion  is  confined  to  the  horizontal  plane,  the  applicable  equations  are  the 
sway  and  yaw  equations.  These  equations  are  standard  Principles  of  Naval 
Architecture,  or  PNA,  equations.  They  are  as  follows; 

Sway: 


Yaw; 


m[v  *ur  ^wp  *f)  -y^ip^  -p)] 

=  (  W  -  B  )  COS0  cos4)  + 


(la) 


m[x^{v  +  ur  +w)  -y^iu  -vr  +  wq)] 

=  ( xw  -  JCpB )  COS0  cos<))  +  ( W  -  ypB )  sin0  +N^ 

If  we  make  the  above  assumptions  in  terms  of  the  variables  in  the  equations: 


1.  No  vertical  motion;  w=0 

2.  W=B 

3.  Small  angle  motion  (linearized  simulation  only) 

4.  No  pitch  rate;  also  negligible  roll;  0=0,  0=0,  (|)=0 

5.  Forward  speed  equals  nominal  speed;  u=U 

6.  Center  of  Buoyancy=Center  of  gravity;  Xg,  yg=Xp,  yp 
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Then  the  sway  and  yaw  equations  become: 


miv  ^  ur  +  xr)  =  Y^ 


1  r  +mx  (v  +  ur)  =  N, 

i  g  / 


The  hydrostatic  forces  on  the  right  hand  side  of  the  equations  are  then  expanded 
in  a  standard  Taylor  series: 


y  =Yz;  +  yr +  Y.t)  +  yf +  y  , 

F  V  r  V  r  totiiitne 


N=Nv  +  Nr+  N.v  +  Nr  +  N  . 

tv  7  V  r  towime 


''^here  y,  =  j  ,  etc, 


Note  that  the  rudder  terms,  and  Yj  have  been  omitted,  as  the  towed  vessel's 
rudder  is  amidships;  there  is  no  applied  force  due  to  the  rudder.  The  only  applied 
force  is  due  to  the  towline.  Substituting  (3a)  into  (2a)  and  (3b)  into  (2b)  and 
rearranging  the  terms  yields: 

Sway: 


(m  -  Y)v-*-i-Y.)r  =  Yv  +  {Y  -mu  )r  +  Y, 


Yaw: 


(I  -  N.  )f  -  N.v  =  N  V  +  N  r  +  N,  .. 

z  r  V  V  T  townte 
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which  can  then  be  rearranged  into  working  form  as 


y  z;  +  ( Y  -m)v  +  (Y^  -mu)\f  +  Y.xif  =  Y  .  (5a) 

r  V  <l>  ~  V~  towline 

Nn+Nu  +  N  W  +  -/  )\i/  =  N  ,  (5b) 

V  V  V  z  '  towltne 

where  r=\}/,  r=Ti)/.  All  that  remains  is  to  describe  the  towline  forces  in  greater 
detail. 


A.  TOWLINE  FORCE 

In  Figure  1,  P  is  the  attachment  point  of  the  tow.  Xp  and  yp  denote  the 
distance  from  the  towed  vessel's  center  of  gravity,  again  using  the  standard  PNA 
body-fixed  coordinate  system.  From  the  diagram,  we  can  see  that  the  horizontal 
force  associated  with  the  towline  in  the  sway  direction  is: 

T  sin  ( Y  +  \|/ ) 

The  moment  causing  the  vessel  to  yaw  is: 

-TXp  ( sin  (y  +  v ) )  "  cos  ( y  +  xy ) 

Assuming  y  and  v[/  are  small  angles: 


sin  ( y  +  \|/ )  =  y  +  ty  = 


y 

L 


=  V 


L 
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Finally,  substituting  for  v,  in  equations  (5a)  and  (5b)  and  rearranging  terms; 
Sway: 


( Y.-m)y  -(Y.u-Y^)\\f 

Ty. 


-  [  Y  w  +  T 


.  i-/ 


(6a) 


Yaw: 


-N^)V 

iv 


-[Nu^Tx 

®  p 


/  \ 

V 


(6b) 


-  —  x  y  *Ty 

p:ip  jp 


These  two  equations  are  motion  describe  the  towed  vessel's  motion.  Note  that  the 
coordinate  y  describes  the  lateral  distance  from  the  towing  path;  yp  describes  the 
distance  of  the  towline  attachment  point  from  the  centerline  of  the  towing  vessel. 
It  is  for  yp  that  the  control  law  is  being  developed. 
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B.  ARRANGEMENT  INTO  STATE-SPACE 

Let  Xi=y,  X2=Vj/,  X3=y,  then,  substituting  into  equations  (6a)  and  6b) 


(Y  -m)x,  +  Y.x. 

V  3  V  4 


=  -Lx+[yu+T 
L  ’  " 


X 

1 

^  L 


(7a) 


N  x.^iN  -I  )x, 

t?  3  y  z  4 


—  XX,  +  [  N  M  +  Tx 

^  P  1  O  P 


\ 

X 

V  ^  J 


X,  -  N  X3  +  ( N^«  -  )  X,  +  _  x^p  +  T 


From  the  definition  of  the  state  variables,  it  can  be  seen  that 


(7b) 


^1=^3 


(8a) 


and 


X  -  X 
•^2  *^4 


(8b) 


These  comprise  the  first  two  of  four  state  equations.  The  second  two  state 
equations  must  come  from  equations  (7a)  and  (7b)  above.  Note  that  these 
equations  are  coupled  in  Xg  and  X4.  Algebraic  elimination  yields  the  final  two 
equations  of  motion.  Multiplying  equation  (7a)  by  (  )  and  equation  (7b) 

by  (  )  yields  an  equation  for  Xg. 
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[(N. -/J(Y, -m)-N,y.]i. 


(  \ 

“ 

Yu+T 

X 

1 

V 

L 

- 

“ 

f  \ 

■ 

Y. 

N  u  +  Tx 

X 

1 

V 

r  p 

L 

- 

- 

*1N  y,-V.(N,-/,)lx3 


+  T 


(N.  -I  ) 
^  2 


f 


(8c) 


Similarly,  multiplying  equation  (7a)  by  (  )  and  equation  (7b)  by  (  m  ) 

and  adding  yields  an  equation  for  . 
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{{Y  -m)  Nu^Tx  1  +_1 
o  ®  p  L 


-N  y  M  +  T  1  +  _I  }r 

1 '  I  i-  JJ  ’ 

.[N.Y  -N,(r-m)]x3 
+  [  -N.  wu  +N,  w  -N,  y.  +N.  V, 


X  N. 

T  (Y  -m)  1+^  -_!  y 
"  L  L  P 


Equations  (8a),  (8b),  (8c)  and  (8d)  can  also  be  expressed  in  state-space  matrix  form 


0  0  1  0 

0  0  0  1  X, 


a  a  0  a  v 

31  32  33  **34  Xj 

a  a  a  a  V 

41  42  43  “44 


-i,,  y. 


Where, 


D  =  (  -  I J  (Y^-m)  -N,Y^ 


+  rx„  1  +  _£ 

^  L 
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a,,=D--  [  (%-IJ  ] 

a3,  =  D--  [  {Y,u-y^}  -Y^(N,u-N^)  ] 

a,.  =D-'-  ^  [  X,  (  y,-m)  -NJ 

1j 


“ 

f 

\ 

1 

“ 

N  ,u  +  Tx^ 

1  + 

- 

Yu  +  r 

1  + 

V  p 

L 

V 

V 

L 

- 

\ 

) 

-* 

- 

\  y-J 

a,,=D---  [  N,y,-N,  {  Y,-m)  ] 

a,,  =D--  [  N,Y^  -  N^Y,  ^  (N^-N,u)  m] 


T 

’  _  V 

1  + 

L 

1  ^  JJ 

=D' 


(y.-m) 


1  + 


N,.. 
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III.  CONTROL  SYSTEM  DESIGN 


A.  SYSTEM  AUGMENTATION 

Equation  (9)  defines  the  state-space  form  of  the  equations  of  motion.  One 
more  state  must  be  defined  and  added  to  this  equation;  that  of  the  ordered,  or 
commanded,  yp.  This  variable  shall  be  called  yp^.  The  usual  response  of  yp  with 
respect  to  time  should  look  like  a  first  order  system,  as  in  Figure  2: 


A 


time 
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The  equation  governing  this  response  is  as  follows: 


T  Y 

p  p 


y  =  Y 

Jp  pc 


-1  1 

or  y  =  —  y  +  —  V 

Jp  'f  'T  ^pc 

P  P 

where  Tp  is  the  non-dimensional  time  constant  for  the  towline  attachment  point 
motion  control  system.  If  yp^  is  now  considered  the  input  to  the  system,  and  yp 
is  considered  a  state  variable,  the  matrix  equation  becomes: 


yp 

-1 

T 

o 

0 

0 

0 

0 

Vp 

'  1  ' 

T 

r 

p 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

^2 

0 

^3 

^3, 

^32 

^33 

®34 

^3 

0 

_ 1 

^41 

^42 

^43 

^44 

0  _ 

Vpc 


(10) 


where  aj,,  ajj,  a^^,  a^,  a^^,  a^,  bj  and  b4  are  as  previously  defined.  In  a 

sense,  the  actual  distance  of  the  towline  attachment  point  yp  is  treated  as  an  extra 
state  to  the  system.  The  control  system  will  be  designed  based  on  this  augmented 
system  of  equations.  The  control  system  will  be  designed  based  on  two  principles 
of  state-space  theory;  Linear  (Quadratic  Regiolator,  and  Sliding  Mode  Control. 


B.  LINEAR  QUADRATIC  REGULATOR 

The  linear  quadratic  regulator  function  arose  from  efforts  to  find  an 
optimum  means  of  control  [Ref.  6].  Normally,  a  gain  matrix  for  a  particular 
control  system  is  calculated  to  achieve  specific  closed-loop  pole  locations.  With 
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the  linear  quadratic  regulator,  or  LQR,  a  specific  performance  criterion  J  (or  "cost 
function")  is  defined,  with  the  only  stipulations  on  the  poles  being  that  they  be 
negative,  or  stable.  This  criterion  J  is  defined  as: 

J  =  j[x'it)Qit)xit)*u'(t)R(t)uU)]dt 

The  matrices  Q  and  R  are  weighting  matrices.  The  matrix  Q  is  a  state 
weighting  matrix,  and  R  is  a  control  weighting  matrix.  The  gains  for  a  control 
system  can  be  calculated  based  on  the  defined  matrices  Q  and  R.  For  instance,  if 
the  elements  of  the  matrix  Q  are  small  relative  to  R,  the  system  will  tolerate  large 
errors  in  the  final  state  with  very  little  control  effort.  Conversely,  if  Q  is  large 
relative  to  R,  very  small  errors  in  the  state  X  will  result,  but  with  considerable 
control  effort. 

C  SLIDING  MODE  CONTROL 

The  second  aspect  of  state-space  theory  to  be  used  in  controlling  the  towed 
vessel  is  that  of  sliding  mode  control.  Since  the  control  law  is  based  on  a 
linearized  set  of  equations  of  motion,  a  lot  of  uncertainty  in  the  response  exists. 
Also,  some  of  the  parameters  may  vary.  A  good  example  for  the  towed  vessel 
system  is  the  towline  tension  T,  which  will  certainly  vary  with  time.  A  control 
law  needs  to  be  developed  that  will  take  into  account  both  the  uncertainties  in  the 
parameters,  and  any  dynamics  that  have  either  not  been  modeled,  or  that  have 
been  linearized.  A  sliding  mode  control  law  can  ensure  both  stability  and 
robustness  of  the  system,  with  the  emphasis  on  robustness.  The  LQR  gains  will 
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ensure  stability  of  the  system.  Sliding  mode  control  is  ideally  suited  to  a  system 
where  the  response  oscillates  between  set  values,  such  as  the  motion  of  a  towed 
vessel;  it  uses  a  high  speed  "toggling"  control  law  to  drive  the  system  onto  a 
desired  "sliding  plane." 

Sliding  mode  control  takes  the  standard  state-space  system 

i:  =  [A]jc  +  [b]M  ^12) 

and  defines  a  sliding  plane 

0(jc)  +S2JC2  +S3X3  (13) 

and  the  coefficient  s,  is  arbitrary.  Equation  (13  )  can  be  written  as 

s^x-O 


where 


s  ^  =  [s,,S2,S3,s^  ] 

Next,  define  the  Lyapunov  function 

V(jc)=i[e(x)]' 

2 

stability  is  guaranteed,  provided  v  (x)  is  a  negative  definite  function.  Another 
way  to  express  this  is 

V(x)=e(:c)6(x)=-Ti'e(x) 
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therefore, 


6  =  -T|^  sign  0 


Since 


e(j:)  =s^(;(:),6(x)  =s'^(i)  =s'^iAx  +bu ) 
so,  substituting  for  0  in  Equation  (15): 


and  solving  for  u: 


s^iAx  +h<)  =  -'n^sign(0) 


u  =  -is'^b)'^  s'^Ax  -V(s^b)’  signiQ) 
Equation  (16)  is  a  sliding  mode  control  law  for  a  generic  system. 


(15) 


(16) 


D.  THE  CONTROL  LAW 

All  that  remains  is  to  define  a  control  law  for  the  augmented  equations  of 
motion  for  the  towed  vessel.  The  control  law  is  defined  using  sliding  mode 
control  and  linear  quadratic  regulator  principles.  Take  the  augmented  system: 
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(17) 


and  split  the  augmented  matrix  A  into  A,„  A,2,  Aj,,  and  A22  as  shown.  Also  split 
the  state  vector  into 


Rewrite  equation  (17)  as 

yi=Aiyi"A2y2"^» 

y2=^2iyi^^22y2"^“ 

The  original  augmented  state  vector  x  now  equals  (y,,  y-^. 
Next,  define  the  sliding  plane  for  the  towed  vessel  as 

e(ar)  =e(y,,yj)=y,  +s’^y2=0 

and 

6  =  -q^sig«(e)  =  y,  +s^yj 


(18) 


(19) 
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but 


and 

so,  substituting  for  y-  and  y-  is  equation  (19): 

+s’"(/ij^y^  *A^y^)  =  -Vi^sign(Q) 

solving  for  u: 

M  =-ll(A„  *s^A^^)y^  -^q^sign(0)  (20) 

This  is  the  system  control  law  for  the  towed  vessel. 

The  gains  [s]  for  the  sliding  mode  control  law  from  minimization  of  the 
linear  quadratic  regulator  cost  function  J.  Recall  that 

/  =  /lUl'[Q)Ul*lul'lRHu))*  <«* 

in  Equation  (21)  the  state  weighting  matrix  Q  is  as  follows: 

0  0  0  0 
0  0  0 

0  0  0  0 
0  0  0  0 

where  qjj  represents  an  error  in  radians  from  the  path  of  the  towing  vessel.  We 
are  choosing  to  weigh  the  state  variable  X2  =  \j/.  The  control  weighting  matrix  R 
will  be  defined  by 
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R  =  (a)-2 

where  a,  the  maximum  non-dimensional  distance  will  be  weighted.  With  Q 
and  R  defined,  the  cost  matrix  function  becomes: 

/  =  min  *Ru^)dt  <22) 

or 

/  =  min  *Ryl)dt  <23) 

The  smaller  the  value  for  R  is,  the  smaller  the  control  effort  required,  but  a  larger 
state  error  in  will  be  have  to  be  accepted  as  a  trade-off.  Conversely,  the  larger 
the  value  for  R  is  the  more  control  required,  with  very  little  state  error  in  q^. 

The  resulting  gains  are  then  placed  into  the  control  law  for  the  augmented 
system  of  equations.  Equation  (20),  and  the  response  of  the  system  is  obtained. 
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IV.  RESULTS 


A.  MATRIXX  AND  FORTRAN  IMPLEMENTATION 

MatrixX  is  an  outstanding  tool  for  understanding  the  response  of  a  control 
system.  A  Fortran  program  was  written  to  use  with  MatrixX.  Since  the  coefficients 
of  the  [A]  and  [B]  matrices  of  the  state-space  equation  are  constant  with  time, 
they  can  be  put  directly  into  the  MatrixX  software,  and  the  resulting  Linear 
Quadratic  Regulator  gains  are  easily  calculated.  Without  the  use  of  MatrixX, 
solving  the  linear  quadratic  integral  is  a  very  difficult  problem  in  numerical 
analysis.  The  actual  simulation  of  the  system  response  was  accomplished  with  a 
Fortran  program  on  the  VAX  computer.  The  LQR  gains  obtained  from  MatrixX 
are  put  into  a  data  file,  along  with  the  non-dimensional  parameters  of  Tp,  LL,  Xp, 
T,  and  the  initial  conditions  of  Y,  \|/,  and  the  maximum  distance  Yp  can  travel 
athwartships.  The  simulation  then  reads  this  data  file,  computes  the  sliding  mode 
control  and  plots  the  response  of  the  parameters  Y,  Yp,  tj/,  and  Yp^  with  respect  to 
time. 

B.  INITIAL  CONDITIONS  AND  CONSTANT  PARAMETERS 

The  results  are  presented  graphically,  and  are  in  non-dimensional  form,  with 
the  exception  of  the  yaw  angle  y,  which  is  in  degrees.  The  non-dimensionalization 
is  standard  Principles  of  Naval  Architecture:  the  relationships  describing  the 
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nondimensional  parameters  are  included  in  the  Appendix.  Two  vessels  were 
studied: 

1.  528  foot  mariner 

2.  1066.3  foot  tanker 

A  summary  of  the  hull  particulars  and  hydrodynamic  derivatives  for  both  these 
vessels  in  also  included  in  the  Appendix.  Four  knots  was  used  as  the  nominal 
forward  velocity.  The  mariner  is  stable  at  four  knots;  the  tanker  is  unstable.  The 
towline  tension  for  both  vessels  is  taken  from  resistance  curves  at  four  knots,  and 
then  non-dimensionalized;  the  value  used  is  0.001  Xp,  the  longitudinal  distance 
from  the  towed  vessel  center  of  gravity  to  the  towline  attachment  point,  is 
assumed  to  be  constant.  A  value  of  0.5  is  used.  The  time  required  for  Yp  to 
"match"  Yp^  as  in  Figure  2,  Tp  is  set  as  0.5.  The  maximum  non-dimensional 
distance  Yp  can  travel  from  port  to  starboard  extremes  is  taken  to  be  0.1;  therefore 
the  maximum  distance  from  an  extreme  to  centerline  is  0.05.  The  initial  conditions 
for  each  data  run  are  Y=0.5,  \|/=0,  and  Yp=0.05  as  in  Figure  3.  Also,  the  weighting 
matrices  were  set  at  q22=5°  (or  .087  radians)  for  (2/  the  state  error  weighting 
matrix,  and  a=0.015  for  R,  the  control  effect  weighting  for  Yp,..  Four  plots  were 
generated  for  each  simulation.  The  y-distance  from  the  path  (Y),  the  towline 
attachment  point  lateral  offset  (Yp),  the  towline  attachment  point  offset  command 
(Yp,,)  and  the  yaw  angle  (y)  are  all  plotted  vs  nondimensional  time.  Variance  of 
parameters  of  interest  are  shown  on  subsequent  nms. 
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1.  Figures  4  and  5:  Standard  Mariner  and  Tanker 

Figure  4  shows  the  response  of  the  mariner  with  the  "standard" 
parameters  and  initial  conditions  stated  above;  Figure  5  shows  the  "standard" 
tanker.  Overshoot  in  Y  and  v  exists  for  both  cases,  although  it  is  higher  in  the 
mariner.  The  Yp^  graph  shows  that  significantly  more  effort  is  required  to  control 
the  tanker,  and  the  oscillation  in  Yp  is  much  greater  than  for  the  mariner.  This 
makes  sense,  as  the  tanker  is  inherently  more  unstable.  Both  vessels  converge  to 
a  straightline  path  in  approximately  60  non-dimensional  seconds.  As  the  tanker 
is  more  unstable,  it  will  be  used  to  demonstrate  all  other  variations  in  parameters. 
Henceforth  all  references  to  "normal"  or  "standard"  tanker  response  will  be  to  the 
response  shown  in  Figure  5. 

2.  Figure  6:  Tanker  with  Y,  i  Not  Observable 

Figure  6  shows  the  vessel  response  if  the  state  variables  Y  and  Y  are  not 
observable  and  are  not  used  in  the  control  law.  Both  Y  and  Y  converge  but  the 
time  required  is  significant;  over  100  seconds.  In  this  case  there  is  no  overshoot 
in  Y,  but  considerable  overshoot  in  \|/.  y,  Yp^  and  Yp  all  oscillate  considerably;  but 
in  the  end,  the  system  does  converge. 

3.  Figures  7,  8:  Observable  Tanker,  Variance  of  6(x) 

Recall  from  Equation  (13)  that  the  sliding  plane  0(x)  is  in  part  arbitrarily 
defined,  as  the  gain  is  chosen  by  the  designer.  For  the  standard  nms,  was 
set  to  be  1.0.  Figure  7  shows  the  tanker  response  if  q^  in  0(x)  is  set  to  0.5,  vdth  all 
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Figure  4.  Mariner  Response  with  Standard  Parameters. 
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Figure  6.  Tanker 


ith  Y,  Y  Not  Observable. 
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state  variables  observable.  The  system  still  converges  in  adequate  time;  the 
parameters  Y,  Yp^  and  Yp  behave  as  in  the  standard  tanker.  There  is 
approximately  twice  as  much  overshoot  in  however.  In  Figure  8,  is  set  to 
2.0;  this  brought  the  overshoot  in  \|f  down  to  the  level  of  the  standard  observable 
tanker,  without  altering  the  favorable  response  of  the  other  three  parameters. 

4.  Figures  9, 10, 11:  Non-Observable  Tanker,  Variance  of  0(x) 

The  state  variables  Y  and  Y  were  assumed  to  be  not  observable  in 
Figures  9  and  10.  In  Figure  9,  r\^  was  set  to  0.5,  with  disastrous  results.  The  towed 
vessel  becomes  unstable.  All  four  parameters  diverge  with  time.  In  Figure  10, 
however,  in  0(x)  is  set  to  2.0,  and  the  system  response  stabilizes.  The  overshoot 
in  V  is  approximately  as  the  same  as  the  standard  tanker,  but  the  overshoot  in  Yp 
and  Yp^  is  twice  that  of  the  normal  tanker.  The  settling  time  in  Y  is  double  that 
of  the  normal  tanker,  but  the  settling  time  in,  Yp^  and  Yp  is  cut  to  forty  seconds. 
Obviously  there  is  a  value  between  T|^  =  0.5  and  =  1.0  in  0(x)  where  the 
response  turns  unstable  for  the  non-observable  case.  In  subsequent  simulations, 
is  set  to  be  2.0.  Figure  11  summarizes  the  response  of  Y  and  Yp  for  the 
observable  and  non-observable,  with  =  2.  The  top  two  graphs  are  for  the 
observable  case,  and  the  bottom  two  describe  the  non-observable  case. 

5.  Figure  13:  Variance  of  the  Sliding  Plane  Switching 

Since  sliding  mode  is  a  high-speed  switching,  or  toggling  method  or 
control,  the  speed  at  which  the  switching  is  to  be  accomplished  must  be  defined. 
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Figure  8.  Observable  Tankers,  ii*-2.0. 
0 
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Figure  9.  Non-Observable  Tanker,  "n^sO-S. 
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Figure  10.  Non-Observable  Tanker,  ti*=2, 


Y  vs  Time 


Figure  11.  Observable,  Non-Observable  Comparison 
of  Y,  Yp  for  t|*=2.0. 


This  switching  speed  can  be  defined  in  terms  of  an  angle.  Note  in  Figure  12  that 
as  the  angle  ^  approaches  90  degrees,  the  system  is  being  asked  to  change 
instantaneously.  Another  way  of  saying  this  is  that  the  system  is  ordering  the 
towline  attachment  point  to  move  instantaneously  from  port  to  starboard.  <|)=90 
degrees  is  impossible;  the  designer  has  to  allow  the  system  time  to  react. 

Figure  13  shows  the  observable  and  non-observable  tanker  response  of 
Y,  Yp  for  (j)=30, 4»=45,  and  <t)=60.  For  the  observable  case,  overshoot,  the  amount  of 
oscillation,  and  settling  time  increase  as  ^  increases.  For  the  non-observable 
tanker,  the  response  destabilizes  between  <j)=45  and  4)=60.  Clearly  there  is  a  limit 
to  the  switching  speed  of  the  system.  (>=45  was  used  for  the  standard  tanker. 

6.  Figiure  14:  Variance  of  LQR  Weighting 

Figure  14  demonstrates  the  effect  of  varying  the  control  effort  in  the 
linear  quadratic  regulator.  As  a  is  increased  from  0.005  to  0.035,  more  effort  is 
used  to  control  the  towed  vessel  motion.  In  Figure  14,  Y  and  Yp  responses  are 
plotted  for  four  different  values  of  a.  In  all  four  graphs,  the  tanker  is  observable. 
Note  that  as  a  increases,  the  overshoot  in  Yp  decreases,  and  the  settling  time 
remains  about  the  same  as  for  the  standard  tanker.  This  makes  sense;  more  effort 
is  being  expended  to  control  Yp  through  Yp^.  The  response  of  Y  remains  essentially 
unchanged. 
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Figure  12.  The  Sliding  Plane. 
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Y.Yp  Y.Yp 


Figure  14.  Tanker  Response  for  Y,  Yp  as  a  Varies. 
_ Yp, _ Y. 
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7.  Figure  15;  Observable  Tanker,  Variance  of  Tp 

Figure  15  demonstrates  the  effect  of  varying  the  time  allotted  for  Yp  to 
match  Ypj  for  the  standard  tanker.  The  Tp  used  for  the  standard  case  was  0.5. 
Varying  Tp  does  not  significantly  affect  the  behavior  unless  the  change  is  in  order 
of  magnitude.  This  result  has  favorable  implications;  if  changing  Tp  doesn't 
radically  affect  the  towed  vessel  resp)onse,  a  smaller  motor  can  be  used  to  drive 
the  device  moving  the  towline  attachment  point,  Yp. 

8.  Figure  16;  Variance  of  Towline  Length 

Figure  16  illustrates  how  changing  the  towline  length  will  change  the 
response  of  the  tanker,  both  in  the  observable  and  non-observable  case.  For  this 
simulation,  the  non-dimensional  towline  length,  LL,  was  lengthened  and 
shortened  by  twenty  percent.  If  the  towline  is  shortened,  the  overshoot  and 
oscillation  in  Y  and  Yp  will  increase  for  both  observable  and  non-observable  cases; 
the  reverse  is  true  is  the  towline  is  lengthened.  Settling  time  will  remain  the  same. 

9.  Figure  17;  Variance  of  Towline  Tension 

Figure  17  shows  how  changing  the  non-dimensional  towline  tension 
will  change  tanker  response.  If  the  tension  in  decreased,  the  overshoot  remains 
about  the  same;  oscillation  and  settling  time  in  Y  and  Yp  increase  for  both  the 
observable  and  non-observable  cases.  The  reverse  is  true  as  the  tension  is 
increased.  A  20  percent  change  in  towline  does  not  cause  a  radical  change  in 
system  response. 
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Figure  16.  Observable,  Non-Observable  Response 
for  Variance  of  LL. 

_ LL=2.0, _ LL=2.5, _ LL=3.0. 
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10.  Non-Linear  (Large  Angle)  Motion 

Recall  from  page  6  that  the  horizontal  towline  force  associated  in  the 
sway  direction  is 

T  sin  ( Y  +  \j/ ) 

and  the  moment  due  to  the  towline  causing  the  vessel  to  yaw  is; 

-TXp  ( sin  ( Y  +  V ) )  -  TypCos  ( Y  +  \|/ ) 

Previously  the  assumption  was  made  that  the  range  of  motion  for  y  and  \|/  was 
less  than  30  degrees.  This  simulation  lifts  this  restriction;  the  sine  and  cosine 
terms  from  the  towline  tension  term  are  left  in  when  running  the  simulation. 
Figure  18  plots  the  standard  tanker  response  for  both  linear  and  non-linear 
tensions.  Note  that  the  tanker  response  remains  stable;  the  sliding  mode  control 
law  is  robust  enough  to  handle  non-linearities  in  the  towline  tensions.  The  non¬ 
line  Fortran  simulation  program  is  included  in  the  Appendix. 
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Figure  18.  Non-Li 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  thesis  explored  the  possibility  of  athwartships  movement  of  the  towline 
attachment  point  as  a  means  of  improving  the  stability  of  towed  vessels.  Newell 
demonstrated  that  use  of  a  full  order  observer  will  work  only  when  the  towed 
vessel  trajectory  is  close  to  that  of  the  towing  vessel  [Ref.  5].  Sliding  mode  control, 
in  conjunction  with  the  Linear  Quadratic  Regulator  cost  function,  provides  a 
much  more  robust  method  of  controlling  the  towed  vessel  motion.  The  sliding 
mode  control  law  can  ensure  stability  of  the  system,  even  if  not  all  the  state 
variables  can  be  observed.  The  only  restrictions  for  the  numerical  applications 
considered  in  this  work  are  that  a  toggling  speed  such  that  ^  is  less  than  an  angle 
between  45  and  60  degrees  is  used,  and  that  the  toggling  gain  ti^  is  greater  than 
a  number  between  0.5  and  1.0.  The  LQR  function  allows  the  designer  to  choose 
where  his  emphasis  is  placed.  By  choosing  whether  to  place  the  emphasis  on  the 
control  effort,  the  designer  can  minimize  the  size  of  the  motor  required  to  drive 
the  mechanical  apparatus  used  to  move  the  towline  attachment  point.  If  the 
designer  wants  to  minimize  steady  state  error,  he  can  change  the  LQR  function 
to  reflect  that  emphasis.  The  LQR  function  can  also  be  used  to  minimize  specific 
state  variables,  such  as  in  the  case  studied,  where  \|/  was  chosen  to  be  minimized. 
It  should  be  emphasized  that  a  "worst-case"  scenario  was  used  to  demonstrate  the 
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robustness  of  the  system;  a  large,  unwieldy  tanker  was  used,  and  two  of  the  four 
state  variables  were  declared  "not  observable."  Even  in  the  worst  case  the  sliding 
mode  control  law  performed  adequately.  Towline  tension  and  length  were 
increased  and  decreased  by  20  percent,  and  the  system  remained  stable.  The  small 
angle  motion  assumption  was  dropped,  and  yet  the  system  continued  to  perform 
satisfactorily. 

B.  RECOMMENDATIONS 

The  sliding  mode  control  law  should  be  robust  enough  to  cover  most  of  the 
unmodeled  dynamics  and  non-linearities,  but  two  of  the  assumptions  made  at  the 
beginning  of  the  thesis  must  be  investigated  further;  the  massless,  inextensible 
towline,  and  the  no  wind,  wave,  or  current  (disturbances)  assumptions.  In  reality 
the  towline  tension  does  not  remain  constant;  consequently,  the  surge  velocity  u 
of  the  towed  vessel  will  also  vary.  A  method  to  take  these  aspects  of  motion  into 
account  needs  to  be  looked  at;  integral  control  is  one  possibility.  Another  is  to 
determine  a  wave  and  towline  function,  and  add  the  two  functions  to  the  right 
side  of  the  sway  and  yaw  equations  of  motion. 

This  thesis  has  shown  that  it  is  possible  to  stabilize  and  control  the  motion 
of  a  towed  vessel  with  a  sliding  mode  control  law  that  uses  linear  quadratic 
regulator  gains.  Further  analysis  needs  to  be  done  on  the  effects  of  disturbances 
and  a  time-variant  towline  tension  on  the  system. 
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APPENDIX 


TABLE  1.  TOWED  VESSEL  DATA 

Vessel 


Property 

Barge 

Taiiker 

Mariner 

LBP,  ft 

191.56 

1066.3 

528 

x'c 

-0.00136 

-0.0009 

-0.0004444 

m' 

0.170 

0.0181 

0.00888 

i'' 

-0.01383 

-0.0171 

-0.00912 

y' 

-0.0153 

-0.0261 

-0.01434 

Jy 

0 

0 

0 

y; 

0.00238 

0.00365 

0.00456 

Ni 

0 

0 

0 

Nv 

-0.007285 

-0.0105 

-0.0046 

H 

1 

0.00188 

0.00222 

0.00115 

-0.00128 

-0.0048 

-0.00296 
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TABLE  2.  NONDIMENSIONAL  TERMS 


m'  =  m  /  irt^ 
u'  =  u  /  U- 
v'  ■=  V  / 

u'  -  u  /  (  ul  -  L) 

v'  ~v  I  (ul-L) 

i'  =  i  /{  ul  -  L)  (  time  ) 

l'  =  I, /  (  m,L^  ) 

r'  =\);'  =\j//  (uj  L) 

r'  -Mf'  =  \if  /  (  u".  /  L‘  ) 

Xg=X^/  1 
Yg^Yg/  L 

Xl  =X.J  (  /  L'  ) 

x',,  =  X,  /  ;n^ 


=  y,  /  (m,u,/L) 

y.'  =  y*  =  y^/ 

y;  =  y^  =  i'y  /  ( ) 
nI^NJ  (m,uj 
/  (  m^L  ) 

Ni  =N^  =W^/  (  ) 

Ni  =Ny  =N^/  (  m.,L^  ) 
T'  =  T  /  {m^ul  /  L) 

1'  =  1/  L 
Xp=Xp/  L 

Y'p^Yp/  L 


(  Tension  ) 
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MATRIX.  LQR  PROGRAM 


Xudot=-0.0009; 

M=0.0181; 

Yvdot=-0.0171; 

Yv=-0.0261; 

Yrdot=0.0; 

Yr=0.00365; 

Ydel=0.0278; 

Nvdot=0.0; 

Nv=-0.0105; 

NrdoUz=-.(X)222 
Nr= -0.0048 
Ndel=-0.0]39 
inquire  LL 
inquire  Xp 
inquire  T 
U=1.0; 

L=1.0; 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

al3=1.0; 

a24=1.0; 

a31=(l/D)*(T/LL)*(NrdotIz-Xp*Yrdot); 

a32=(l/Dr((Yv*U+T*(l-t-Xp/LL))*NrdotIz-Yrdot*(Nv*U-(-T^Xp*(l+Xp/LL))); 
a33=(l  /D)*(Nv*Yrdot-Yv*NrdotIz); 
a34=(l/D)*((Yvdot*U-Yr)‘NrdotIz-Yrdot*(Nvdot*U-Nr)); 
a41={l/D)*(T/LL)*(Xp*(Yvdot-M)-Nvdot); 

a42=(l/Dr((Yvdot-M)*{Nv*U-t-T*Xp*(l-f-Xp/LL))-Nvdot*{Yv*U+'r(l-t-Xp/LL))); 

a43=(l/D)*(Nvdot*Yv-Nv*(Yvdot-M)); 

a44=(l/D)*(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U)*M); 

b3=(l /D)»T*(NrdotIz/LL-Yrdot»(l +Xp/LL)); 

b4={l/D)*T((Yvdot-M)*(l-hXp/LL)-Nvdot/LL); 

A=[0,0^1 3m0,0^24;a3U32^33,a34;a41  ,a42,a43,a44]; 

B=l0,O;b3;b4l; 

C=I0A1,0]; 

D=[01; 

Q=[0  0  0  0;0  131.3316  0  0,0  0  0  OO  0  0  01; 

inquire  alpha 

R=[(l/alpha)”21; 

<ElG,K>=REGULATOR(A,B,Q,R); 

INQUIRE  TP 

Al=l-l/TP,0,0,0,0;0,0,0,l,0O,0,0,0,l;B3,A31,A32,A33,A34;B4,A41,A42,A43,A44l; 

B1=[1/TP;0,-0,-0;0]; 

S=(1,K1; 

G1=-INV(S*B1)*S»A1; 

G2=-INV(S*B1); 

ETA2=G2 

GAIN=G1 

SIGMA=S 
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LINEAR  SIMULATION  PROGRAM 


REAL  KYP,Kl,K2,K3,K4,L,LL,Nvdot,M,Nv>lr,Nrdot 

Xudot=-0.0009 

M =0.01 81 

Yvdot=-0.0171 

Yv=-0.0261 

Yrdot=0.0 

Yr=0.00365 

Ydel=0.0278 

Nvdot=0.0 

Nv=-0.0105 

Nrdot=-0.00222 

Nr=-0.0048 

Ndel=-0.0139 

U=1.0 

L=1.0 

OPEN  (10,nLE='LINEAR.DAT',STATUS='OLD') 
OPEN  (ll,FILE='Y.DAr,STATUS='NEW') 

OPEN  (12,F1LE='PSI.DAT',STATUS='NEW') 

OPEN  (13,nLE='YP.DAT',STATUS='NEW') 

OPEN  (14,nLE='YPC.DAT',STATUS='NEW') 

READ  (10,*)  TSIM,DELTAT,IPRNT 
READ  (10,*)  KYP,K1,K2,K3,K4 
READ  (10/)  SYP,SLS2,S3,S4 
READ  (10/)  GEAR 
READ  (10/)  TP,LL,XP,T 
READ  (10/)  ETA2,PH1 
READ  (10/)  Y,PSI,SAT 

1SIM=TSIM/DELTAT 

PI=4.0*ATAN(1.0) 

PHI=PHITI/ 180.0 

PSI=PSI*PI/180.0 

SATP=  SAT 

SATM=-SAT 

XI  =Y 

X2=PS1 

X3=0.0 

X4=0.0 

YP=0.0 

YPC=0.0 


D=Nrdot*(Yvdot-M)-Nvdot*Yrdot 

al3=1.0 


nnnn  non  n  non 


a24=1.0 

a31  =(1  /D)‘(T/LL)*{Nrdot-Xp*Yrdot) 

a32=(l /Dr((Yv*U+'r(l +Xp/LL)rNrdot-Yrdot*(Nv*U+'rXp*(l +Xp/LL))) 
a33=(l  /D)*(Nv*Yrdot-Yv*Nrdot) 
a34=(l/D)'^((Yvdot*U-Yr)*Nrdot-Yrdot*{Nvdot*U-Nr)) 
a41=(l/D)*(T/LL)*(Xp*(Yvdot-M)-Nvdot) 
a42=(l /D)*((Yvdot-M)*(Nv*U+T*Xp*(l +Xp/LL)) 

&  -NydofCYv^U+Td+Xp/LL))) 
a43={l  /D)*(Nvdof Yv-Nv*(Yvdot-M)) 
a44=(l/D)*(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U)*M) 
b3=(l/D)*T»(Nrdot/LL-Yrdot*(l+Xp/LL)) 
b4=(l/D)*T*((Yvdot-Mr(l+Xp/LL)-Nvdot/LL) 

SIMULATION  BEGINS 

DO  1  I=1,ISIM 
YPDOT  =-YP/TP+YPC/TP 
XIDOT  =  X3 
X2DOT  =  X4 

X3DOT  =  B3*YP  +  A31*X1  +  A32*X2  +  A33*X3  +  A34*X4 
X4DOT  =  B4*YP  +  A41*X1  +  A42*X2  +  A43»X3  +  A44*X4 

YP  =  YP  +  YPDOTDELTAT 
XI  =  XI  +  XI  DOr  DELTA! 

X2  =  X2  +  X2IX)T*DELTAT 
X3  =  X3  +  X3DO'PDELTAT 
X4  =  X4  +  X4DOrDELTAT 

CONTROL  LAW 

SIGMA  =  SYP’YP  +  SrXl  +  S2*X2  +  S3*X3  +S4»X4 
SATSGN=  SIGMA/TAN(PHI) 

IF  (SATSGN.GE.{+].0))  SATSGN=+1.0 

IF  (SATSGN.LE.(d.O))  SATSGN=-1.0 

UHAT  =  KYP‘YP  +  KrXl  +  K2»X2  +  K3*X3  +  K4*X4 

UBAR  =  ETA2*GBAR»SATSGN 

YPC  =  UHAT+UBAR 

IF  (YPC.GE.SATP)  YPC=SATP 

IF  (YPC.LE.SATM)  YPC=SATM 

Y  =X1 
PSI=X2 


PRINT  RESULTS 


I=J+1 

IF  (J.NE.IPRNT)  GO  TO  1 
TIME=I*DELTAT 
WRITE  (*  •)  TIME,Y 
WRITE  (11,*)  T1ME,Y 
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WRITE  {12,»)  TIME,PS1*1 80.0/PI 
WRITE  (13,*)  TIME,YP 
WRITE  (14,*)  TIME,YPC 
1=0 

1  CONTINUE 
STOP 
END 


NON-LINEAR  SIMULATION  PROGRAM 


C  NONLINEAR  SIMULATION  -  TANKER 
C  SURGE  NOT  INCLUDED 
C 

REAL  KYP,K1,K2,K3,K4,L,LL.NVD0T,NV,IZ,NRD0T,NR,NVRR 

c 

XUDOT=-0.0009 
M  =  0.0181 
YVDOT=-0.0171 
YV  =-0.0261 
YRDOT=  0.0 
YR  =0.00365 
NVDOT=  0.0 
NV  =-0.0105 
IZ  =0.0 
NRDOT=-0.0022 
NR  =-0.0048 
YVRR  =-0.045 
NVRR  =  0.0061 
UTOW  =  10.0 
L  =  1066.3 
RHO  =  1.9905 
SB  =680625.0 
P  =  9.78 
Q  =  1.93 

SB  =  SB/(0.5*RHO*L*L*UTOW*UTOW) 

C 

OPEN  (10,nLE='LINEAR.DAT',STATUS='OLD') 

OPEN  (ll,nLE='Y.RES',STATUS='NEW') 

OPEN  (12,nLE='PSI.RES',STATUS='NEW') 

OPEN  (13,nLE=TP.RES',STATUS='NEW') 

OPEN  (14,nLE='YPC.RES',STATUS='NEW') 
OPEN(15,nLE='NONLINEARl.RES'^TATUS='NEW') 

C 

READ  (10,*)  TSIM,DELTAT,IPRNT 
READ  (10,*)  KYP,K1,K2,K3,K4 
READ  (10,*)  SYP,S1,S2,S3,S4 
READ  (10,*)  GBAR 
READ  (10,*)  TP,LL,XP,T 
READ  (10,*)  ETA2,PHI 
READ  (10,*)  Y,PSI,SAT 
C 

ISIM  =  TSIM/DELTAT 
PI  =  4.0*ATAN(1.0) 

PHI  =  PHI*PI/180.0 
PSI  =  PSPPI/180.0 
SATP  =  SAT 
SATM  =-SAT 
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non  n  n  n  n  non  non 


X  =  LL-XP 

u  =  1.0 

V  =  0.0 
R  =0.0 
YP  =0.0 
YPC  =0.0 

c 

DEN=(M-YVDOT)*(lZ-NRDOT)-YRDOT*NVDOT 
SIMULATION  BEGINS 
DO  1  I=LISIM 

F2=YVV+0.5*YVRR*V*R»R+{YR-MPR 
F3=NV*V+0.5*NVRR*V*R*R+NR*R 

SING=(Y+XP*SIN(PSI)+YP*COS(PSI))/LL 
GAMMA=ASIN{SING) 

YPIXDT  =-YP/TP+YPC/TP 
VDOT  =(F2-rSIN(GAMMA+PSI))/(M-YVDOT) 

ROOT  =(F3-'DXP*SIN(GAMMA+PSI)-T*YP*COS(GAMMA+PSI))/(IZ-NRDOT) 
YDOT  =U*SIN(PSI)+VCOS(PSI) 

PSIDOT=R 

YP  =YP  +  YPDOT  •DELTA! 

U  =U  +UDOT  ’DELTA! 

V  =V  +VDOT  ’DELTA! 

R  =R  +RDOT  ’DELTA! 

X  =X  +XDOT  ’DELTA! 

Y  =Y  +YDOT  ’DELTA! 

PSI  =  PSI  +  PSIDOTDELTAT 

CONTROL  LAW 

SIGMA  =  SYP*YP  +  Sl’Y  +  S2*PSI  +  S3’YDOT  +  S4*PSIDOT 
SATSGN=  SIGMA/TAN(PHI) 

IF  (SATSGN.GE.(+1.0))  SATSGN=+1.0 
IF  (SATSGN.LE.(-1.0))  SATSGN=-1.0 

UHAT  =  KYP*YP  +  Kl’Y  +  K2’PSI  +  K3*YDOT  +  K4’PSIDar 
UBAR  =  ETA2’GBAR’SATSGN 
YPC  =UHAT+UBAR 
IF  (YPC.GE.SATP)  YPC=SATP 
IF  (YPC.LE.SATM)  YPC=SATM 

PRINT  RESULTS 


J=I+1 

IF  (J.NE.IPRNT)  GO  TO  1 
TIME=I*DELTAT 
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WRITE  (V)  TIMELY 

WRITE  (ll,*)  TIME,Y 

WRITE  (12,*)  TIME,PSI*180.0/PI 

WRITE  (13,*)  TIME,YP 

WRITE  (14,*)  TIME,YPC 

WRITE  (15,*)  TIME,Y,PSI*180.0/PI,YP,YPC 

1=0 

1  CONTINUE 
STOP 
END 
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